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If dark matter is composed of neutralinos, the gamma-ray radiation produced in their annihilation offers an 
attractive possibility for dark matter detection. This process may contribute significantly to the extragalactic 
gamma-ray background (EGB) radiation, which is being measured by the FERMI satellite with unprecedented 
sensitivity. Using the high-resolution Millennium-II simulation of cosmic structure formation we have produced 
the first full-sky maps of the expected contribution of dark matter annihilation to the EGB radiation. Our 
maps include a proper normalization of the signal according to a specific supersymmetric model based on 
minimal supergravity. The new simulated maps allow a study of the angular power spectrum of the gamma-ray 
background from dark matter annihilation, which has distinctive features associated with the nature of the 
annihilation process. Our results are in broad agreement with analytic models for the gamma-ray background, 
but they also include higher-order correlations not readily accessible in analytic calculations and, in addition, 
provide detailed spectral information for each pixel. In particular, we find that color maps combining different 
energies can reveal the cosmic large-scale structure at low and intermediate redshifts. 
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1. Introduction 

Although dark matter accounts for most of the mat- 
ter in the Universe, its nature remains unknown and 
its presence has only been inferred through its grav- 
itational effects. However, if dark matter is made of 
neutralinos, a new particle predicted by Supersym- 
metry, it would also interact (although very weakly) 
with ordinary matter, and it might be detected soon 
in laboratories on Earth. 

In addition, neutralinos can self-a nnihilate to pro- 
duce ordinary par ticles like positro ns Baltz and Edsiol 
p"999j ■ neutrinos iBerezinskv et alj jl996{ and gamma- 
ray photons. If these byproducts of the annihilation 
are copious enough, they could be detected soon. In 
the present work, we will focus on the gamma-ray pho- 
tons as a residual of the annihilation process. 

This gamma-ray radiation is expected to be pro- 
duced more abundantly in regions with high dark 
matter density. Thus, it seems best to look for 
it in very dense nearb y regions, such as the cen - 



ter of our own Ga l axy Bertone and MerrittI [2005l | 



iJacholko wska et "aL |2006{" and/or th e centers of its 
satellite galaxies Wood et al.l |2008| . Actually, it 



turns out that the best prospects for the detection of 
gamma-rays from our Galaxy are obtained by looking 
slightly off-center to avoid confusion of the signal with 
other sources of gamma r a ys residing at the Ga lactic 
center IStoehr et al.l |2003| , ISpringel et al.l [2008 1. 

Outside our Galactic halo, gamma-rays are also pro- 
duced in large quantities by the annihilation of dark 
matter in all the many halos and subhalos within 
our past light-cone, contributing to the so-called ex- 
tragalactic gamma-ray background (EGB) radiation. 
The EGB has been measured by different satellites, 
in particular in the energy range be tween 1 MeV 
and 30 GeV by COMPTEL and EGRETlStron£ et al.l 
(2004| . Our understanding of the EGB is greatly im- 



which covers an energy range between 20 MeV and 
300 GeV and features an improved sensitivity com- 
pared with its predecessor EGRET. 

Although the EGB also rec eives co ntributions from 
other sources, such as blazars lAndo et al. [2007a] and 
cosmic rays accelerat ed at structure formation shocks 
lJubelgas et al.l [2008j |. the energy spectrum and an- 
gular power spectrum of the annihilation radiation 
have distinctive features that m ay open up effective 
ways fo r disentangling the signal Aiido and Komatsul 
' 2006l| . ISiegal-Gaskins and Pavlidoul |2009f . There- 



proving by the Large Area Tele scope aboard the re - 
cently launched FERMI satellite lAtwood et al.l|2009j |. 



fore, a detailed analysis of the EGB is a viable possi- 
bility for detecting dark matter. 

Previous studies have analyzed this possibility us- 
ing analytic approaches. However, it is not clear 
how accurately these methods capture the non-linear 
structures resolved in the newest generation of high- 
resolution cosmological simulations. Also, the previ- 
ous analysis have so far been restricted to statistical 
statements about the power spectrum of the EGB, or 
its mean flux. For a full characterization of the sig- 
nal, it would be useful to have accurate realizations 
of maps of the expected gamma-ray emission over the 
whole sky. 

In our study I Zavala et al.l |2009j |. we therefore 
focus on predicting the EGB radiation directly 
from cosmological N-body simulations. To this 
end, we use the Mi l lenniu m-II simulation (MS-II) 
iBovlan-Kolchin et al.l [20091 . With 10^" particles in 
a homogeneously sampled volume of (100 /i~^Mpc)"^, 
it is one of the best resolved structure formation simu- 
lations to date. We use a ray-tracing technique to ac- 
cumulate the signal from all halos and subhalos over 
the past light-cone of a fiducial observer, positioned 
at a plausible location of the Milky Way in the simu- 
lation box. 

The values we calculate for each pixel of such a sky 
map are that of the specific intensity, the energy of 
photons received per unit area, time, solid angle and 
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energy range: 



A,o = ^ j £7.0(^^7,0(1 + z),z)^ 



dr 



(1) 



where the integral is over the whole line of sight, r 
is the comoving distance and E^^ is the energy mea- 
sured by the observer at z = 0. The gamma ray emis- 
sivity (energy emitted per unit energy range, unit 
volume and unit time) associated with dark matter 
annihilation is given by: 



{av) 



dE^ 



Px 



(2) 



where dN^/dEj is the total differential photon spec- 
trum summed over all channels of annihilation, 
and are the mass and density of neutralinos, and 
{av) is the thermally averaged product of the annihi- 
lation cross section and the M0ller velocity. Note that 
e-yfi is evaluated at the blueshifted energy (1 -I- z)Ej^o 
along the line-of-sight to compensate for the cosmo- 
logical redshifting. 
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Figure 1: SUSY factor as a function of 7-ray emission energy 
for a selection of benclimark points representatives of tlie 
allowed regions in the mSUGRA parameter space: CA 
(black), FP (red), RAF (green ) and B (blue). 



2. Dark matter annihilation and the SUSY 
factor 

Of all the quantities in Eq. ([2]), only depends 
on the spatial distribution of neutralinos, the rest is 
related to its intrinsic properties, which we can con- 
veniently analyze separately by defining the so called 
SUSY factor: 



/sUSY = 



dN-y {av) 
dE^ TO^ 



(3) 



In the present section we describe the particular SUSY 
model that we used to compute the gamma-ray emis- 
sivity produced by dark matter annihilation. 

We restrict our analysis to the framework of the 
minimal supergravity model (mSUGRA), which is 
popular, among other reasons, due to its relative sim- 
plicity. The large number of free parameters in general 
SUSY is reduced to effectively four free parameters in 
mSUGRA: mi/2, mo, Aq, which are the values of the 
gaugino and scalar masses and the trilinear coupling, 
all specified at the GUT scale, tan/3 which is the ratio 
of the expectation values in the vacuum of the two 
neutral SUSY Higgs, and finally, the sign of fi, the 
Higgsino mass parameter. 

The general 5-dimensional parameter space of the 
mSUGRA model is significantly constrained by var- 
ious requirements: consistency with radiative elec- 
troweak symmetry breaking (EWSB) and experimen- 
tal constraints on the low energy region. It is also nor- 
mally assumed that if the major component of dark 
matter is the lightest neutralino, then its relic den- 
sity fl^ should be equal to the observed abundance of 



dark matter rioM- This condition reduces the allowed 
regions in the parameter space considerably. 

The different regions in the mSUGRA parameter 
space that satisfy these constraints have been studied 
often in the past and have received generic names: (B) 
bulk region (low values of toq and TO1/2), (FP) focus 
point region (large values of toq), (CA) co-annihilation 
region (low mo and m-^ ^ mf^ , where ti is the lightest 
slepton) and (RAF) rapid annihilation funnel region 
(for large values of tan/3 and a specific condition for 
m^^). Instead of making a full scan of the allowed pa- 
rameter space, we define several "benchmarks points" 
which fulfill all the experimental constraints discussed 
above and which can be taken as representative of 
the different regions that we just described. The 
benchmar k points were sel ecte d following the analy- 
sis of Battagha et al.l (2004j and lGondolo et al.l l2004l 
We u s e the numerical code DarkSUSY iGondolo et al.l 
j2004 , 12005 ^1 to analyze these benchmark points. 

Fig. [T] shows the energy spectrum of /susy com- 
puted for four of the benchmark points analyzed, 
which are representatives of the different regions: 
CA (black), FP (red), RAF (green ) and B (blue). 
The main features of Fig. [T] are determined by the 
photon spectrum dN^/dE^ which receives contribu- 
tions from three mecha nisms of photon production 
iBringmann et all 2008| : (i) gamma-ray continuum 
emission following the decay of neutral pions produced 



^The relation which should hold is: 2m^_Ri rn^ where A is 
the CP-odd Higgs boson; see for example lFentJ |2005ll for details. 
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during the hadronization of the primary annihilation 
products; (ii) monoenergetic gamma-ray hnes for neu- 
trahno annihilation in two-body final states contain- 
ing photons; (iii) internal bremsstrahlung (IB), which 
leads to the emission of an additional photon in the 
final state; the IB contribution to the spectrum is typ- 
ically dominant at high energies resulting in a charac- 
teristic bump. Processes (ii) and (iii) are subdominant 
to process (i), but they display distinctive spectral fea- 
tures intrinsic to the phenomenon of annihilation. 

In Fig.m we have shown only 4 of all the benchmark 
points we have considered; the rest have very similar 
features and lie in between these 4 cases. Fig. [T] also 
indicates that the normalization of the spectrum of 
/suSY typically varies by three orders of magnitude 
among the different regions of the allowed parameter 
space. 

For the rest of our analysis, we will use one par- 
ticular benchmark point, which is the one we high- 
light with a thick blue line in Fig. [TJ It is a model 
with — 185 GeV that gives an upper limit on 
/suSY for the benchmark points in the bulk region 
and is also close to the maximum value of /susy for 
all the benchmark points analyzed. The latter is the 
main reason to choose this particular model because 
it is desirable that the predicted gamma-ray flux has 
a large value among the different theoretical possi- 
bilities. There are other reasons as well. The mass 
of the neutralino in this model is "safely" larger than 
the lower mass bounds coming from experim ental con- 
strain ts {m^ > 50 GeV according to Heister et al.l 
20041 1. but low enough to be detectable in the ex- 
periments available in the near future. Also, a high 
value of the parameter tan/3 seems to be favored by 
other theoretical expectations iNunez et al.l [2008j | . 



3. The dark matter annihilation 
contribution to the EGB 

The specific intensity (Eq. [Ij describes the total 
emission from dark matter annihilations integrated 
over the full backwards light cone along a certain di- 
rection. We use the data of the MS-II to fill the whole 
volume contained in the past-light cone of an observer 
located at a fiducial position at z — 0. Since there are 
only 68 simulation outputs in total, we can approx- 
imate the temporal evolution of structure growth by 
using at each redshift along the past light-cone the 
output time closest to this epoch. To cover all space, 
we use periodic replication of the simulation box. 

Our strategy is not to account for the emission at 
the level of individual dark matter particles, but in- 
stead to use entire halos and subhalos, allowing us 
to accurately correct for resolution effects. The list 
of dark matter particles at each output time is hence 
replaced by a catalog of dark matter substructures. 



We further assume that these substructures are well 
represent ed by a sp herically symmetric NFW density 
profile .Navarro et al, .1997.] . Under this assumption 
the gamma-ray luminosity fro m a given substructur e 
can be written as a scaling law lSpringel et aL 1 |2008j : 

/I 23 
p'NFwir)dV^-—^, (4) 
^ 'max 

where Vinax is the maximum rotational velocity of the 
subhalo and r^^x is the radius where this maximum is 
reached. The most recent high-resolution N-body sim- 
ulations support a slightly revised density profile that 
becomes gradually shallower towards the center, the 
so called Einasto profile. For this case, the coefficient 
1.23 in Eq. (|4]) changes to 1.87, increasing our results 
effectively by 50%. This is an insignificant difference 
compared to other uncertainties in the analysis of the 
dark matter annihilation radiation. For the purposes 
of our work, the NFW profile is a good approxima- 
tion and using it simplifies comparisons to a number 
of results in the literature. 

Our map-making method then becomes a task to 
accurately accumulate the properly redshifted emis- 
sion from these structures in discretized representa- 
tions of the sky and to include a resolution correc- 
tion for poorly resolved objects. All our maps use 
^pix = 12 X 512^ ^ TT X 10^ pixels, corresponding to 
an angular resolution of ^ 0.115°. 

A given pixel in the simulated maps covers a solid 
angle Afipix. Our map-making code computes the 
average value of the specific intensity within the area 
subtended by this solid angle by conservatively dis- 
tributing the emission of each substructure over the 
appropriate pixels. Combining Eqs. ([2]-|4l) we can 
write the relevant sum over all substructures in the 
light-cone as: 

^ ^ Lhw{dh,rh)EjfifsijsY{zh)\E.y,oi (5) 

where the function w^dh, r^) is a weight function that 
distributes the luminosity of a given halo onto the pix- 
els overlapping with the projected "size" of the halo; 
the latter depends on the distance of the observer to 
the halo dh and the transverse distance rh between 
the halo center and the center of the pixels it touches. 
Except for structures that are very nearby, the high 
central concentration of the emission of a subhalo and 
the limited angular resolution of our maps give most 
subhalos the character of unresolved point sources. 

For the purpose of extending the predictions of our 
maps down to the damping scale limit of neutrali- 
nos 1O~^M0), we divide the maps into separate 
components: (i) the contribution of resolved halos and 
subhalos with a minimum mass of 6.89 x 10^ /i~^Mq 
(called "ReHS" in the following); (ii) the contribu- 
tion of unresolved main halos with masses in the range 
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1.0 X 10-6-6.89 xIO^/j-IMq (referred to as "UnH"); 
and finally, (iii) the contribution of unresolved subha- 
los in the same mass range as in the case (ii) (compo- 
nent "UnS"). 



3.1. Resolved halos and subhalos 
(ReHS) 

In Fig. [5J we show full sky maps of the the 7-ray 
emission of all main halos and subhalos that are de- 
tected in the MS-II; this corresponds to an essentially 
perfectly complete sample of all halos above a mass 
limit of 6.89 X 10^ /i^^Mq. In the top panel of Fig. H a 
partial map at low redshift (corresponding to the first 
shell at z = Q) is shown in order to illustrate typical 
foreground structure, whereas the bottom panel gives 
an integrated map out to z = 10, which is approxi- 
mately the full EGB from annihilation. The maps in 
Fig. [2] are for _E-y.o — 10 GeV. In the color scale used 
for the maps, red corresponds to the highest and black 
to the lowest values of the specific intensity. 

3.2. Unresolved halos (UnH) 

The maps shown in Fig. [2] are complete down to 
the minimum mass in the MS-II that we can trust, 
^ 6.89 X 10^ /i~^Mo. In order to make a prediction of 
the full EGB coming from dark matter annihilations, 
we extrapolate the 7-ray flux to account for the con- 
tribution of all missing dark matter halos down to the 
cutoff mass 10"^ /i~^M0. For this purpose we use ex- 
trapolations of the power law behavior found for the 
total gamma-ray luminosity of host halos as a func- 
tion of mass. This is an analysis that wb present in 
detail in section 4.1 of I Zavala et al.l 2009j |. 

The way we incorporate this extrapolated contribu- 
tion in the 7-ray maps is the following. We assume 
that the EGB radiation from the missing halos in the 
mass range 10"^ /i^^M© to 6.89 x 10** /i^^M© is 
distributed on the sky in the same way as the one 
from the smallest masses we can resolve in the sim- 
ulation, which we adopt as the mass range between 
1.4 x 10^ /j-iM© and ~ 6.89 x 10^ /i-iM©. Hence, we 
compute the value of a boost factor bh with which each 
resolved halo in the mass range 1.4 — 6.89 x 10^ /i~^Mq 
needs to be multiplied such that the luminosity of the 
unresolved main halos is accounted for as well. We 
found that bh ~ 60, and that its value is nearly inde- 
pendent of redshift up to the highest redshift we can 
reliable make the extrapolation, which is z ~ 2.1. 



3.3. Unresolved subhalos (UnS) 

To add the contribution of unresolved substructures 
to the 7-ray maps, we first assume that the identified 
subhalo population is complete down to a mass limit 



of Msuh ~ 6.89 X 10^ h^^yiQ. Below this mass we 
compute the gamma-ray luminosity of subhalos using 
extrapolations of the power law behavior of the con- 
tribution of these subhalos to the luminosities of their 
host halos. The paramet ers of the power law are given 
in detail in section 4.2 of IZavala et al ] (2009f . We add 
this contribution in different ways for the following 
two cases. 

In the case of resolved main halos with more than 
100 particles, we distributed the extra luminosity com- 
ing from unresolved substructures among the subhalos 
of the main host, assuming in this way that unresolved 
subhalos are distributed in the same way as the re- 
solved ones. If the main halos have no subhalos, the 
extra luminosity is given directly to them. 

For the second case of all main halos with masses 
less than 6.89 x 10* /i^^Mq, we use the following strat- 
egy. From the results described in the last subsection 
we can get the total luminosity coming from all main 
halos with masses between the damping scale limit 
and 6.89 x 10* /i~^M0. Thus, we can compute the 
boost factor &sub for which this luminosity needs to 
be multiplied in order to include the full contribution 
of all subhalos of the main halos in this mass range. 
We found that 6sub is roughly independent of redshift. 
Using the range of values found for the power laws de- 
scribed before, we obtain that bsuh ~ 2 — 60. We take 
the extreme values of this range in the following results 
as the minimum and maximum values that reflect the 
uncertainties in our extrapolation method and should 
bracket the true result. 



3.4. Isotropic and anisotropic 
components of the EGB from dark 
matter annihilation 

The upper panel of Fig. [3] shows the mean value of 
the specific intensity A/^^o divided by the comoving 
thickness of the shell at £'^^0 — 10 GeV as a function of 
redshift. Recall that each partial map represents the 
total specific intensity in a shell of constant comoving 
thickness, so this provides information about where 
the signal is coming from. The black line is for the 
ReHS case, while the dashed line is for main halos only 
(resolved and unresolved, including the boost factor 
bh = 60). The remaining lines are for the resolved and 
unresolved main halos and their subhalos, boosted to 
include the contribution of structures all the way down 
to the damping scale limit; here the blue line is for the 
"minimum boost" with b^ub = 2 whereas the red line 
is for the "maximum boost" with &sub = 60, but both 
include the boost bh — 60 for unresolved main halos. 

In all cases, the contribution from partial maps up 
until z ^ 1 is almost constant, indicating that the 
increasing number of sources for shells at higher red- 
shifts, due to the larger volume seen behind each pixel, 
is approximately compensated by the distance factor 
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Figure 2: Upper panel: A partial map showing the EGB produced by dark matter annihilation in nearby structures. Only 
sources within 68 Mpc from the observer are considered for the map. The color scale gives a visual impression of the values of the 
specific gamma-ray intensity for each pixel in the map; the red color corresponds to the highest values of specific intensity. The 
observed energy of the simulated gamma-ray radiation is 10 GeV. Lower panel: The full gamma-ray sky map from dark matter 
annihilation containing sources up to 2 ^ 10. In both maps only the contribution of the smallest halos resolved by the simulation 
(~ lO" h~^'!sAQ) has been taken into account. 



and the spectral effects from the cosmological redshift- 
ing and the intrinsic variation of the emission spec- 
trum, as described by /susy (see Fig. [1]). 

The lower panel of Fig. [3] shows the accumulated 
mean value of the specific intensity as a function of 
redshift. This clearly shows that the most relevant 
contributions come from maps up to z 2. Overall, 
applying the maximum boost increases the mean value 
of of the ReHS case by three orders of magnitude. 
For most energies, the /susy spectrum is monotoni- 
cally decreasing. In these regions, the dependence on 
Ij,o with redshift will have the generic form shown in 
Fig. |3l However, at the highest energies the shape 
is expected to change dramatically due to the impor- 
tance of Internal Bremsstrahlung and/or of monochro- 
matic lines (see Fig. [IJ close to the rest mass energy 



of the dark matter particle. Such an effect is clearly 
visible in the black dotted line in Fig. [31 which is the 
ReHS case for an energy of 100 GeV and features a 
bump in the redshift distribution, a reflection of the 
importance of IB at the highest energies. 

The angular power spectrum of the EGB is an im- 
portant tool to study the statistical properties of its 
anisotropy on the sky. In fact, certain classes of 7- 
ray sources are expected to exhibit different power 
spectra, making this a potential means to identify the 
origin of the unresolved EGB. 

In Fig. 21 we show the angular power spectrum 
l{l + l)C//27r as a function of the multipole at an 
observed energy of E^^ — 10 GeV, for all the cases 
discussed in Fig. [31 At large scales, I ^10, the power 
spectrum is related to the clustering of dark matter 
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Figure 3: Upper panel: Mean annihilation intensity per 
comoving siiell tiiickness, as a function of redshift for 
individual partial maps. The solid black line is for resolved 
halos and subhalos. The dashed black line is for main halos 
only, resolved and unresolved down to the cutoff mass using a 
boost factor bf^ = 60. The blue and red solid lines show the 
contribution of all components, resolved and unresolved halos 
and subhalos, boosted with the extreme values in the interval 
fesub = 2 — 60. All lines are for -£^7,0 = lOGeV, except for the 
black dotted line which is for the ReHS component with 
^^7,0 = 100 GeV. Lower panel: The same as the upper panel 
but for the accumulated intensity, /7,o(> 1 + 2). 



halos. When only the main halo contribution to the 
EGB is considered (dashed black line), the normaliza- 
tion is lower because most of the 7-ray signal comes 
from low mass halos that are less clustered than more 
massive halos. 

The blue line in Fig. |4l corresponding to the full ex- 
trapolation including subhalos but with the minimal 
boost 6sub = 2, has exactly the same power spectrum 
than the case with main halos only, at all scales. This 
is because the signal from main halos is dominant in 
this case, and subhalos have a negligible effect. In 
contrast, in the case where subhalos have a signifi- 
cant contribution mediated by 5sub — 60 (red line), 
the normalization is larger at small scales because the 
signal is dominated by subhalos belonging to the most 
massive halos, and the latter are strongly clustered. 

As the angular scale decreases, I > 10, the power 
spectrum depends more and more on the internal 
structures of halos. For the cases where substructures 
are ignored or are negligible (dashed-black and blue 
lines), the slope becomes steeper, with a slope close to 
2. The power spectrum for the cases where substruc- 




10-^1 ^^^^^ ^^^^^ ^^^^^ 

1 10 100 1000 

multipole I 



Figure 4: Angular power spectrum of the EGB produced by 
dark matter annihilation as a function of the multipole / for 
£'7,0 = 10 GeV. The line styles are as in Fie. [51 except for the 
dotted lines that show the predictions from lAndo et aO 
|2007bl for a subhalo-dominant contribution with and without 
considering tidal destruction (lower and upper dotted lines), 
respectively. 

tures are relevant (red and black solid lines) behaves 
differently, however. In the range I £ [20, 100], it be- 
comes slightly shallower, i.e. the signal is slightly more 
isotropic in this regime. Contrary to the strong cen- 
tral concentration of the matter in a halo, the number 
density profile of subhalos is considerably shallower 
than a NFW profile and produces a lu minosity profile 
in pro jection which is essentially flat ISpringel et al.l 
[2008]. This effect continues until I - 200 where the 
power spectrum becomes dominated by the low-mass 
main halos. 

The dotted line s in Fig. [H were taken from the an- 
alytical results of lAndo et a l. 2007b] and are shown 
here f or comparison, see t he las paragraphs of section 
5.5 in I Zavala et al.l j2009j ] for a discussio n on the dif- 
ference s between our results and those of lAndo et al.l 
|2007bl| . 

3.5. Energy dependence 

An important result that we obtain from the anal- 
ysis of our results is that there are differences in the 
shape of the power spectra for different energies. We 
can exploit these differences and enhance the power 
at small or large angular scales by taking the ratio of 
sky maps at different energies. The upper panel of 
Fig. [5] shows the total sky maps for E^^q = 0.1 GeV 
and Ej Q — 32 GeV in the left and right, respectively, 
smoothed with a Gaussian beam with a FWHM of 
5°. Both maps show only small anisotropics, with the 
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Figure 5: Upper panel: Full-sky maps at energies 0.1 GeV and 32 GeV in the left and right, respectively. The maps were 
smoothed with a Gaussian beam with a FWHM of 5°. At a single energy, a full-sky map is very smooth, nearby structures are 
only minimally visible. Lower panel: Ratio of the maps in the upper panel (left) and a partial map containing only nearby 
structures within 68 Mpc for an observed energy of 0.1 GeV (right). Creating difference maps ("color" maps) using different 
energy channels greatly enhances the signal of nearby structures. 



exception of a couple of prominent structures. The 
sky map in the lower left panel is the ratio of the 
two maps; it clearly enhances the signal of nearby 
structures. Convincing supporting evidence for this 
enhancement comes from the nearest partial map at 
z = for i^^^o = 0.1 GeV shown in the lower right 
corner of Fig. [SJ Most of the prominent dark matter 
structures that can be seen in this map are also clearly 
present in the map with the energy ratio. We conclude 
that the 7-ray sky maps that the FERMI satellite will 
obtain at different energies could be used to construct 
difference maps which may then show enhanced cor- 
relations with nearby cosmic large-scale structure. 
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